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ABSTRACT 

We study the impact of lens galaxy environments on the image separation distribution 
of lensed quasars. We account for both environmental convergence and shear, using a 
joint distribution derived from galaxy formation models calibrated by galaxy-galaxy 
lensing data and number counts of massive elliptical galaxies. We find that the ex- 
ternal field enhances lensing probabilities, particularly at large image separations; the 
increase is ~30% at 9 = 3" and ~200% at 9 — 5", when we adopt a power-law source 
luminosity function $>(L) oc L . The enhancement is mainly driven by convergence, 
which boosts both the image separation and magnification bias (for a fixed lens galaxy 
mass). These effects have been neglected in previous studies of lens statistics. Turning 
the problem around, we derive the posterior convergence and shear distributions and 
point out that they are strong functions of image separation; lens systems with larger 
image separations are more likely to lie in dense environments. 
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1 INTRODUCTION 

Strong gravitational lensing offers a unique probe of cos- 
mology and the physical properties of early-type galaxies. 
For example, the ensemble properties of a well-defined sam- 
ple of stron g lens systems can constrain the cosm ologi- 
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cal constant lTurnerlll990l: iFukueita et alj|l99cl; iKochaneld 
Il996bl : IChaell2003l : iMitchell et al]l2005lh the density profile 
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and evolution of early-type galaxies IIKochanek et al. 2000; 
iRusin et all l2003allbt iRusin fc Kochanekl 120051). and the 
amount of substructure in lens galaxies ijMgtcalf fcM adaul 
l200ll:lDalal fc Kochanekll2002t Kochanek fc DalaJl2004) . In 
addi tion, the image separation distribution of lensed quasars 
jTurner et al.lll984) probes the efficiency of baryon cooling 
inside dark matter halos of different masses, and thereby 
constrains models of galaxy form ation llKochanek fc White! 
l200ll: lKeetonll2001al: lQgurill2002l) . So far, more than 80 lens 
systems have been discovered, and the largest homogeneous 
subsample suitable for statistical studies is the Cosmic Lens 
All-Sky Survey (CLASS), containing 13 l ensed radio sources 
at image separatio ns 0'.'3 < 6 < 3" iMvers et alJ l2003t 
iBrowne et al.l20 03h A larger statistical sample is being con- 
structed from the Sloan Digital Sky Survey (SDSS), which 
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has already discovered more than 10 new lenses i ncluding the 
largest-separation lense d quasar known to date lllnada et alJ 
l2003UQguri et alJl2004) . 



Lensing probability distributions are often computed 
using a clean spherical lens object without external fields. 
However, in reality, most observed lenses have non-spherical 
galaxies, and many require a significant external tidal shear 
iKeeton et al Jl997l) . While a portion of the shear ma y come 
from fluctuations of matter along the line o f sight iSeliakl 
Il994t lBar-KaneJll996t iMomcheva et al]l2005l) . much of it is 
thought to be associated with mass in the immediate en- 
vironment of the lens galaxy. Most lens galaxies are early- 
type, w hich lie prefer e ntially in dens e environments. For in - 
stance, iKeeton et alJ i2000t see also iBlandford et al.ll200ll) 
used galaxy demographics (specifically, the galaxy lumi- 
nosity function as a function of type and environment) to 
predict that at least ~25% of lens galaxies lie in groups 
or clusters of galaxies. Indeed, more than 20 lens galax- 
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Figure 1. Top: Image separation distributions of all lensed 
quasars (solid) and those whose lens galaxies are confirmed to 
or suspected to lie in a group or cluster (dashed). Note that not 
all lensed quasars have had their environments characterized ob- 
servationally. Bottom: The fraction of lensed quasars in groups 
and clusters as a function of image separation. 
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Conventional wisdom holds that ellipticity and en- 
vironment mainly affect the relative numbers of dou- 
ble and quadruple lenses and do not significantly 
modify the total lensing probability <|Ko chanek| |l996b: 
iKeeton et al.lll997tlRusin fe Tegmarkll200lHChadl200l but 
se^ lKeeTonfcZal3ludoflll2o7)3r^ showed 
explicitly that ellipticity and shear change the lensing proba- 
bility by only a few percent for most sour ce luminosity func- 
tions (see also iPremadi fc Martell I2004T) . They also found 
that shear (and to a lesser extent ellipticity) shifts and 
broadens the distribution of image separations for a given 
lens galaxy. However, they did not put the two pieces to- 
gether and discus s changes in th e overall image separation 
distribution. Also. lHuterer et alJ (and everyone else, for that 
matter) neglected the effects of external convergence from 
the environment. External convergence is often omitted from 
models of individual lenses be cause the mass sheet d egen- 
eracy renders it unmeasurable iGorenstein et aDll988l : rSahal 
2000). Nevertheless, theoretical studies show that external 
convergence affects lensing analyses in many important ways 
iKeeton fc Zabludoflll20oil . 

There is already some observational evidence that lens 
galax y environments af fect the image separation distribu- 
tion, ioguri et alJ i2005l) found an overabundance of lensed 
quasars with image separations 9 ~ 3". They argued that 
the excess arises from lens galaxy environments, since many 
of the 9 ~ 3" lenses appear to lie in groups or clusters. 
We extend this consideration to all lensed quasars: Fig- 
ure shows the image separation distribution for all lensed 
quasars and for the subset whose lens galaxies lie in groups 
or clusters. Lens galaxy environment is clearly correlated 



with image separation such that lenses with larger separa- 
tions tend to lie in groups and clusters. The correlation is al- 
most tautological at large separations (say 9 > 4"), because 
it is hard to achieve a large deflection angle without some 
environmental convergence. But the figure shows a gradual 
increase of the fraction even at 9 < 4". This trend provides 
further evidence that environmental converge nce and shear 
are im portant: it contradicts the prediction bv lKeeton et alJ 
(2000) that, if external convergence and shear are irrelevant, 
lenses in dense environments should have a smaller mean 
image separation than lenses in the field (because dense en- 
vironments tend to have a larger ratio of dwarf to giant 
galaxies than the field). 

While these data certainly suggest a connection between 
environment and lensing, the observational evidence is in- 
complete because syste matic surveys of lens e nvironments 
have only just begun fsee lMomcheva et al.l20o'5T) . It is there- 
fore valuable to undertake a theoretical analysis of how en- 
vironments affect lensing probabilities, and what range of 
environments is expected for lenses with various image sep- 
arations. In this paper, we compute the lensing effects of 
shear and conve rgence associated w ith matter near the lens 
galaxy (see also iMoller et alJl2002t for the effect of groups 
on lens statistics). We use a realistic model for the distri- 
bution of galaxy environments, which is based on TV-body 
simulations and the halo occupation distribution and cali- 
brated by observations of galaxy-galaxy l ensing and num- 
ber c ounts of massive elliptical galaxies iDalal fc Watsor] 
2005). For the lens objects, we assume a singular isother- 
mal sphere (SIS) because it is both simple and a reason- 
able model for the density profile of lens galaxies (e.g., 
[Cohn et_jjj|200jyTreu_fcI^c ^m et alJ 

boost iRusin fc Kochanekl2005h . We neglect the ellipticity of 
lens galaxies for simplicity, but note that the effects of el- 
lipticity are generall y smaller than those of external shear 
iHuterer et alJl2005h . 

This paper is organized as follows. Section |2] describes 
the joint probability distribution function of convergence 
and shear that we use throughout the paper. Section |3| re- 
views the calculation of lensing probabilities and the image 
separation distribution. In 2] we examine the image sepa- 
ration distribution from various viewpoints in order to un- 
derstand how lens environments affect the distribution. Sec- 
tion summarizes our main conclusions. Throughout the 
paper we assume a concordance cosmology with mass den- 
sity Om = 0.3, cosmological constant Qa = 0.7, and dimen- 
sionless Hubble parameter h = 0.7. 



2 JOINT PDF OF CONVERGENCE AND 
SHEAR 

We consider the probability distribution functions (PDFs) of 
convergence and shear originating from lens galaxy environ - 
ments, using the model derived bv lDalal fc Watsonl (2005). 
They placed galaxies in an iV-body simulation using a halo 
occupation distribution calibrated to match number counts 
and tangential shear profil es measured for mass ive elliptical 
galaxies in the SDSS (see ISheldon et al] 120041) . and deter- 
mined the convergence and shear at each galaxy position. 
Figure 3 of their paper shows the resulting distributions. The 
convergence K ex t and shear |7ext| have similar means » 0.03, 
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Figure 2. Radial profile of convergence and shear in a typical 
massive dark matter halo. We assume an NFW halo with mass 
M v i r = 10 14 h~ 1 Mq and concentration parameter c = 5.5. The 
lens has redshift Z; = 0.45, while the source has redshift z s = 2. 



but very different tails: large convergences (ftext > 0.1) are 
more common than large shears (|7ext| > 0.1). There are 
two reasons for the difference. First, for a typical group 
or cluster dark matter halo, we expect that a large shear 
is always accompanied by an even larger convergence. Fig- 
ure [5] illustrates this idea by showing the convergence and 
shear profiles fo r a dark matter halo modeled with an NFW 
density profile iNavarro et al.l fl997T) : in the central region 
(k, \ j\ > 0.1), we see that k > |7|. 1 Second, convergence and 
shear sum differently when there are multiple perturbers. 
The probability distribution of convergence is highly skewed, 
with a tail extending to large positive values but no corre- 
sponding tail at large negative values. In contrast, the dis- 
tributions of both shear components are symmetric about 
zero. Thus, a large convergence provided by one perturber 
cannot be canceled by contributions from other perturbers, 
while large shear values can be canceled. As we shall see, the 
longer tail of convergence is one reason that convergence is 
more important than shear. 

Figure|5]also implies that convergence and shear are not 
independent but rather correlated. Therefore in this paper 
we consider the joint PDF of convergence and shear, which 
is shown in Figure|3] The correlation between K ext and |7ext| 
is quite evident, particularly for large values. 

Although the environment distribution was derived for 
a lens redshift zi = 0.45, we assume little evolution in 
the elliptical galaxy population over the range < z < 1 
(see ISchade et alJ Jl99a llm et alJ 120021 : lOfek et all l2003t 
IChae fc Maoll2003ft . Thus, we simply extrapolate the sur- 
face mass density associated with environment to all red- 
shifts. Since convergence and shear depend on the surface 



1 More generally, for a power-law density profile p(r) <x 
the convergence and shear are related as \ r y\/n = \rj — 1|/|3 - 



Therefore the argument holds for different mass profiles, as long 
as they have inner profiles that are shallower than the isothermal 
sphere (rj = 2). 
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Figure 3. Contours of the joint PDF p(/c ext , |7 e xtl) of external 
convergence and shear derived bv lDalal & : Watson (2005), assum- 
ing redshifts z t = 0.45 and z a = 2. Contours are spaced by 0.5 
dex, with thicker lines indicating higher probabilities. 



mass density in units of the critical density for lensing, 
E CI .i t oc D oa / (D iDi s ), the redshift dependence of the PDF 
for convergence and shear is straightforward. The PDF was 
derived for galaxies with a velocity dispersion o ~ 216 km 
s _1 , which is typical for lenses with image separation 6 ~ 1", 
but we apply it to all galaxies. In other words, we neglect any 
correlation between a galaxy's velocity dispersion and envi- 
ronment. However, we know that galaxies with higher ve- 
locity dispersion will tend to reside in denser environments; 
accordingly, our estimate of environmental effects on wide- 
splitting lenses is strictly a lower limit. We exclude extreme 
cases with re oxt + |7cxt| ^ 1, because those would be catego- 
rized as lensing by groups/clusters rather than by a galaxy 
perturbed by its environment. 



3 LENS STATISTICS WITH EXTERNAL 
CONVERGENCE AND SHEAR 

We model each lens object as an SIS galaxy with external 
convergence and shear. The Einstein radius of an SIS with 
velocity dispersion o is 



V C / Doa 



(1) 



where D oa and Di a are, respectively, angular diameter dis- 
tances from the observer and lens to the source. 

The important properties of the lens are the image sep- 
aration as a function of source position, the lensing cross 
section, and the magnification bias. We define the image 
separation 8 to be the maximum distance between multiple 
images, which is a well-defined and observable quantity. It 
is useful to define the normalized separation 6 such that 



0(u) = E 0(u) 



(2) 
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Note that 6(u) = 2 for all u in the absence of convergence 
and shear. 

We combine the cross section and magnification bias to 
compute the "biased cross section" , 



BA = 



*(L) 



1 du, 



(3) 



where the integral is over the multiply-imaged region of the 
source plane, and we assume a power law source luminosity 
function 4>l{L) oc L - ' 3 . Unless otherwise specified, we adopt 
the source luminos ity function of the CLASS survey, j3 = 2.1 
jMvers et all2003h . We take /i to be the total magnification 
of all images, which is appropriate for surveys in which lenses 
are identified from high-resolution follow-up observations of 
unresolved targets (including both CLASS and SDSS). We 
choose coordinates such that the source position u is in units 
of 9e, which means that BA is naturally (and conveniently) 
in units of 6 E . 

We compute image separations an d biased c ross sec tions 
using the public software gravlens bv iKeetonl (l2001bl) . For 
each set of (/text, I Text |) we place 10 random sources in the 
smallest circle enclosing the lensing caustics. We solve the 
lens equation to find the images, and then compute the im- 
age separation and total magnification. Finally, we sum over 
the multiply-imaged sources to compute the biased cross sec- 
tion. 

The next step is to integrate over lens galaxy popula- 
tions to obtain the total lensing probability. This involves 
integrating over appropriate distributions of galaxy masses 
(or velocity dispersions), redshifts, and environments: 



7 cdt . 
dzi——{l + zi) 

dzi 



dKext d 1 7ext | P ( «cxt , 1 Text | ) 



/, dn 
da- (Do, 



9e) BA. 



(4) 



Here D a i is the angular diameter distance from the observer 
to the lens, and p(/t e xt, |7ext|) is the joint PDF of exter- 
nal convergence and shear from §5] We specify the distribu- 
tion of galaxy velocity dispersions using the velocity func- 
tion dn/da of early- type galaxies determ ined from ^3 0,000 
galaxies at 0.01 < z < 0.3 in the SDSS JSheth et alJl2003t 
iMitchell etaIll2005lL 

The image separation distribution of lensed quasars can 
then be obtained by differentiating equation Q: 



dP 
U 



. 2 d(BA) 



d §l£-p(D ol 6v 
Q d6-E da 



I- 



dO 

dKextd|7exi 

a dn (D ol 0y d(BA) 



(•») 



cdt 

dzi——(l + zi 

dzi 



|p(«ext, I Text I) 



29 da e 2 
We fix the source redshift to z 3 
redshift for lensed quasars. 



(6) 



d& 

= 2.0 which is a typical 



4 RESULTS 

4.1 Dependence on convergence and shear 

Before presenting results that account for full distribution of 
convergence and shear, it is useful to study how fixed values 
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Figure 4. Effects of external convergence and shear on the pre- 
dicted image separation distribution. Here we assume certain 
fixed values of convergence and/or shear instead of taking ac- 
count of full joint PDF. Top: We vary the convergence ft e xt from 
to 0.5 (bottom to top), while the shear is fixed at |7cxt| = 0. 
Middle: We vary the shear |7ext| from to 0.5 (bottom to top at 
9 > 1"), while the convergence is fixed at fcappa ex t = 0. Bottom: 
We set K cx t = |7cxt| and change the value from to 0.5 (bottom 
to top at 9 > 1"). 



of Kext and/or |Text| affect the image separation distribution. 
These results are shown in Figure [1] 

Convergence and shear both modify the shape of image 
separation distribution, particularly at relatively large image 
separations. For instance, at 6 — 5" a convergence K cxt = 0.2 
increases the lensing probability by one order of magnitude, 
and a shear of |7 cx t| =0.2 enhances it by a factor 2-3. For 
the same values, convergence clearly has much more effect 
on lensing probabilities than shear. Combined with the fact 
that large values of convergence are more likely to occur than 
large values of shear (see , we expect that convergence has 
the greater impact on image separation distributions. 

The results can be understood as follows. First, exter- 
nal convergence magnifies the image plane with respect to 
the source plane, which increases the image separation by 
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Figure 5. Image separation distributions with (solid) and with- 
out (dashed) environmental convergence and shear, calculated 
from equation We now use the joint distribution of exter- 
nal convergence and shear from Figure 131 
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Figure 6. Enhancement of the lensing probability as a function of 
image separation. The enhancement is defined to be Pext/-Po — 1, 
where P ex t is the probability with the external field while Pq is 
the probability without. We show the enhancement for both the 
differential lensing probability dP/dd and the cumulative proba- 
bility P(> 6). 



a factor (1 — ftcxt) -1 . Second, while external convergence 
does not affect the caustics, it does change the biased cross 
section by giving additional magnification to the images. 
Specifically, the biase d cross section increases b y a factor 
(1 - Kcxt)" 2 '' 3 " 1 ' fsee lKeeton fc Zabludof3l2004 . As a re- 
sult of these two effects, convergence shifts the image sep- 
aration distribution up and to the right in Figure 2] The 
increase is significant where the image separation distribu- 
tion is a deceasing function (9 > 1"), but negligible where it 
is an increasing function (9 <C 1"). By contrast, shear mainly 
broadens the image separation distribution. It modestly en- 
hances the lensing probability at 9 > 1", and suppresses it 
at 9 < 1". 



4.2 Full result 

We are now ready to consider the image separation distribu- 
tion with the full joint PDF of convergence and shear, which 
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Figure 7. Enhancement of the lensing probability, computed at 
an image separation 8 = 3", as a function of the slope j3 of the 
source luminosity function. The vertical line indicates the fiducial 
value we adopt in this paper (f) = 2.1, as appropriate for the 
CLASS lens survey). 



is shown in Figure The enhancement of the lensing prob- 
ability, which is plotted in Figure |S| is small at 9 ~ 1", but 
significant at larger image separations: the increase is ~30% 
at 9 — 3" and ~200% at 9 — 5". Thus, lens galaxy environ- 
ments are very important when we discuss the shape of the 
image separation distribution. Indeed, they can account for 
the ex cess of lensing probabilities reported bv lQeuri et all 

In Figure HJ we also plot the enhancement of the cumu- 
lative lens probability P(> 9). We find that the enhance- 
ment of the total lensing probability integrated over the all 
image separation is ~ 15%. This enhancement is equivalent 
to that obtained by shifting the cosmological constant by 
AQa ~ 0.05, therefore it cannot be ignored in accurate de- 
termination of cosmological parameters with lens statistics. 

As discussed above, some of the effects of convergence 
and shear are mediated by magnification bias: convergence 
directly magnifies images, while shear generates and length- 
ens the tangential caustic and thus increases the cross sec- 
tion for high magnifications. The enhancement of the lensing 
probability is therefore expected to depend on the source 
luminosity function. Indeed, Figure |7| shows that the en- 
hancement at 9 — 3" increases strongly with the luminosity 
function slope (3 (recall <& oc L~°). This result implies that 
bright optical lensed quasars with large image separations 
are more likely to lie in dense environments. 



4.3 What enhances the lensing probability? 

Dissecting the results further helps us understand the 
changes to the image separation in more detail. As discussed 
in H4.ll we believe that convergence is more important than 
shear because of its stronger effect on the lensing probability 
(see Fig. and the longer tail to high ft cxt in the joint PDF 
(see Fig. |HJ. To check this hypothesis, we project the joint 
PDF to a distribution of converge nce or shear alone (the 
same distributions shown in Fig. 3 iDalal fc Watsonl 120051) 
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Figure 8. Top: Image separation distributions that include only 
convergence (dashed), only shear (dotted), or both (solid). Bot- 
tom: Contributions to the image separation distribution from 
different re ox t intervals. From top to bottom (at 9 ~ 1"), the 
dotted curves correspond to K cxt < 0.05, 0.05 < K cx t < 0.15, 
0.15 < Kcxt < 0.35, and 0.35 < fc cx t. The full image separation 
distribution is shown with the solid line for reference. 



and recompute the image separation distribution. The re- 
sults, shown in the upper panel of Figure |5] confirm that 
the enhancement of the lensing probability is in fact driven 
by convergence. 

In that case, it is useful to understand what values of 
K cxt contribute most to the lensing probability. We see this 
by decomposing the image separation distribution into con- 
tributions from different fi: cxt intervals, as shown in the bot- 
tom panel of Figure |H| The interpretation depends on the 
image separation. At small image separations (9 < 3"), most 
lenses are expected to have small convergences (/text J$ 0.05), 
and larger and larger convergences are more and more rare. 
At intermediate separations (3" < 9 < 4"), the largest con- 
vergences (« ex t > 0.35) become more important than in- 
termediate values. Finally at large separations (8 > 5"), 
the largest convergences become dominant. In other words, 
while dense environments with large convergences may be 
rare, nearly all of the largest separation lenses will be found 
there. 



4.4 Environments of lens galaxies 

Figure |S] suggests that the typical environments of lenses 
are very sensitive to image separation. Put another way, 
the distribution of convergence and shear for actual lens 
galaxies will differ from the distribution shown in Figure 
13 (which applies to normal elliptical galaxies), by a factor of 
the lensing probability. To quantify this effect, we derive the 
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Figure 9. Posterior distributions of convergence (top) and shear 
(bottom), taking into account lensing probabilities. Distributions 
for lenses with image separations 8 = 1" (solid) and 3" (dashed) 
ares shown. The original (unweighted) distributions at zi = 0.45 
are shown by dotted lines for reference. 



posterior PDF of convergence and shear for lenses at fixed 
image separations of 9 — 1" or 3". Figure [5] shows that the 
distributions differ from each other, and from the original 
(unweighted) distributions. The distributions for lenses with 
9 = 1" are similar to the original distributions, just shifted 
to slightly larger values. However, the distributions for lenses 
with 9 = 3" show particularly large increases in the tail of 
high convergence or shear values. 

To see this more clearly, in Figure Hol we show the frac- 
tion of lenses with a "strong" environment (defined to be 
K cx t or |7 cx t| > 0.1), as a function of image separation 8. 
The fraction increases monotonically as the image separa- 
tion grows. For instance, at 9 = 1" the fraction of lenses 
with a strong convergence is 11%, while the fraction with a 
strong shear is 6%. By 8 — 3" the fractions have increased 
to 24% and 13%, respectively. (For comparison, the frac- 
tions of non-lens elliptical galaxi es with strong converge nce 
or shear are just 6% and 3%; see lDalal fc Watsonll2005L 'l In 
other words, the distribution of lens galaxy environments is 
a strong function of image separation, and this effect should 
be taken into account when comparing observed environ- 
ments with theoretical predictions. The increasing fraction 
of "strong" environments with image separation is qualita- 
tively consistent with the data shown in Figure 
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Figure 10. Fraction of lenses with "strong" environments (de- 
fined by K C xt or 1 7oxt | > . 1 ) , as a function of the image separation 
9. 

5 SUMMARY AND DISCUSSION 

We have studied how the convergence and shear from lens 
galaxy environments affect lens statistics, in particular the 
distribution of lens image separations. We find that the ex- 
ternal field enhances the lensing probability, especially at 
large image separations, and the effect increases with the 
slope of the source luminosity function. We argue that the 
enhancement is driven mainly by the external convergence, 
which has been neglected in p revious studies of th i s sort . 
Our results mesh with those of iKeeton fc Zabludofj (120041) 
to indicate that it is essential to include convergence from 
the environment in order to obtain correct results from both 
lens modeling and lens statistics. For example, since the ex- 
ternal field changes the shape of the image separation distri- 
bution rather than its overall amplitude, it will bias attempts 
to constrain the velocity dispersion function of early-type 
galaxies us ing the lens image separation distribution (e.g., 
IChadl2005l) . 

The environmental boost in the lensing probability de- 
pends on the image separation, which means that the poste- 
rior distribution of lens environments does as well. The frac- 
tion of lenses with a "strong" environment (a convergence 
or shear larger than 0.1, say) increases monotonically with 
image separation. Large separation lenses are more likely to 
be found in dense environments. This bias will be impor- 
tant when comparing the observed distrib ution of lens en- 
viron ments to theoretical predictions (e.g.. lMomcheva et alJ 
120051) . 

One puzzle remains: even when we take the separa- 
tion/environment correlation into account, we find that the 
predicted fraction of lens galaxies with a large shear |7cxt| > 
0.1 is ~10%. However, models of mo st observed 4-imag e 
lenses require shears of ~0.1 or larger iKeeton et al.lll997l) . 
While 4-image l enses are certainly biased toward dense en- 
vironments (see Holder & Schcchter] | 2003l for a theoretical 
perspective, and [Momch^va^t^dTl200Rl for intriguing obser- 
vational results), it does not appear that the bias is strong 
enough to reconcile our predicted shear distribution with 



the "observed" values. At present, it is not clear whether 
the problem lies with the lens models or the predicted en- 
vironment distributions. If the latter, then our conclusions 
regarding the importance of lens galaxy environments should 
clearly be revisited. 

It is also possible that relaxing some simplifying as- 
sumptions in our analysis will affect the quantitative results. 
The most important simplification is the assumption that 
all galaxies have t he same distribution of convergence and 
shear. As shown in [Sheldon et alJ l)2004T) . however, the tan- 
gential shear signal on scales < 1 Mpc increases with 
increasing velocity dispersion, suggesting that more massive 
ellipticals reside in denser environments than less massive 
ellipticals. We expect this effect to enhance the importance 
of environment for wide-splitting lenses relative to the esti- 
mates presented here. 

Another simplification was our neglect of lens galaxy 
ellipticity. The effe cts of ellipticity are somewhat similar to 
those of shear (see lHuterer et aTT l2005'l. which has a mod- 
est effect on the lensing probability. Therefore we expect 
that our results would not change significantly with the ad- 
dition of ellipticity. Nevertheless, since ellipticity is impor- 
tant in predictions of the relative numbers of 4-image and 
2-image lenses, it will need to be considered more carefully 
when comparing observed and predicted environment distri- 
butions for 4-image lenses alone. 

It seems to be too early to say whether there is con- 
flict between the observed and predicted distribution of lens 
galaxy environments. But it is an intriguing question that 
deserves to be studied both observationally and theoreti- 
cally. 
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